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SUMMARY 


If in radiative equilibrium due to adiabatic expansion a 
rising mass of gas becomes warmer than its surroundings, its 
motion will continue. Under this condition convection currents 
are stable. Therefore, convection currents can be produced only 
if the amount of the adiabatic temperature gradient is smaller 
than the radiative gradient. The convective heat flux can be com- 
puted from the energy balance during convection. It is also pos- 
sible to derive the convective flux from the equation of radiative 
transfer, which applies to non- equilibrium photon transport pro- 
cesses, by substituting convective for the radiative parameters. 


Symbol 

a 
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c v 

E(r) 
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H(r) 
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LIST OF SYMBOLS 
Definition 

Radiation density constant 
Source function 
V elocity of light 

Specific heat at constant pressure per unit mass 

Specific heat at constant volume per unit mass 

Energy density (energy per unit volume) 

Acceleration due to gravity at distance r from 
the center 

Radiative power flux-density at distance r from 
the center 

Radiant power per unit area per unit solid 
angle at r 

Emitted power per unit mass per unit solid 
angle at r 


Symool 

l 

L(r) 

P(r) 

p < r >rad 

r 

R* 

S 

T(r) 

U 

v 

7 

e 

~K 

P(r) 

a 

T 

do> 

V 

Subscripts 

ad 

k 

rad 

Superscript 


LIST OF SYMBOLS (Concluded) 

Definition 
Mixing length 

Radiant power at the distance r from the center 

Hydrostatic pressure at the distance r from 
the center 

Radiation pressure at the distance r from 
the center 

Radius from the center of the Sun to some point 
of the interior 

Universal gas constant 

Scale height 

Temperature at the distance r from the center 
Internal energy per unit of convective mass 
V elocity of a stable convective element 
Ratio of specific heats, Cp/c v 

Angle between the direction of a ray and the axis 
of a small cylinder 

Rosseland's K (corrected for induced emission) 
Mean molecular weight of the Sun's material 
Density at the distance r from the center 
Stefan- Boltzmann constant 
Optical depth 
Element of solid angle 

Logarithmic gradient, temperature to pressure 

adiabatic 

convective 

radiative 

prime: refers to convecting element 
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I. INTRODUCTION 


The physical properties of the Sun are explained by a com- 
posite model for which the energy transport mode is radiative in 
the core and radiative and convective in the envelope. If the con- 
dition of radiative equilibrium is maintained, no unsteady con- 
vective mass motions can occur. However, if the conditions for 
radiative equilibrium are not fulfilled, perturbed mass and energy 
transport may result. The basic conditions which affect stability 
are stated by the four differential equations for hydrostatic equilib- 
rium, conservation of mass, luminosity, and temperature gradient 
(convective and radiative). These equations have to be comple- 
mented by three constitutive equations in order that the behavior 
of the gases may be characterized. The constitutive equations 
are the equation of state, the equation for the absorption coefficient, 
and the equation for energy generation by nuclear processes. 

The temperature gradient has two distinct forms, adiabatic 
and radiative. Which form is applicable is determined by whether 
the stability condition, 


dT 

dr 


< - (1 - 


J_. T_ dP 
y’ P dr ’ 


is met (adiabatic) or is not met (radiative). 

By using only the convective energy transport condition, 
the entire outer region of the Sun is described by the adiabatic 
gradient. This is a rough procedure which yields erroneous re- 
sults. Since both radiative and convective fluxes occur in con- 
vective equilibrium, both contributions are considered in each of 
the two derivations outlined in the next two sections. 


a 
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H. DERIVATION OF THE CONVECTIVE FLUX 
FROM THERMODYNAMIC EQUILIBRIUM 


For the calculation of energy transport by radiation and 
convection in the convective layer of the Sun, the pure adiabatic 
and radiative gradients are supplemented by an overall average 
gradient and a gradient describing a turbulent element in sur- 
roundings which are otherwise in radiative equilibrium. These 
four gradients, negative in the sense of decreasing temperature 
with increasing radius r, are in order of decreasing magnitude: 

V r ad = (dlnT/dlnP) rac j for radiative equilibrium 

V = (dlnT/dlnP) actual average gradient at 

radius r 

V' = (dlnT/dlnP) / for an average convective 

element 

V ad = (dlnT/dlnP) a ^ for adiabatic changes of state. 

If as a result of adiabatic expansion a mass of gas rises 
through surroundings which are otherwise in radiative equilibrium, 
and if the negative adiabatic temperature gradient of the rising 
mass is less than the negative radiative gradient of the surround- 
ings, then the rising mass will cool less than the surroundings 
and thus be relatively warmer. If expansion exactly compensates 
for relative warming, then stable mass motion will continue. 

These conditions, adiabatic expansion with temperature gradient 
of the rising mass less than the radiative gradient of the surround- 
ings, are the necessary and only conditions for generation of stable 
convection currents. The gradient relationship can be expressed 
in the form: 

^ V'<v< v rad . 
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The convective power density can be computed by the 
following procedure. When a turbulent mass rises a distance Ar, 
its temperature drops less rapidly than that of its surroundings 
and thus the mass tends to be warmer than its surroundings: 

AT = [(dT/dr)' - (dT/dr)] Ar. 

(The primed and unprimed terms have the same meaning as in 
the gradient relationships in the preceding paragraph. ) 

The convectively transported power density equals: 

H k = c p p AT v = c p pT(AT/T)v. 

(The bars above the symbols indicate mean values over a layer of 
thickness dr. ) 

The formula for hydrostatic pressure is determined as 

follows: 

dP = -gpdr 

dlnP = ~(g)i /R*T)dr= -dr/S. 

The mean ascent of an arbitrary turbulent element is set equal to 
half the mixing length: 

Ar = dr = , 

and the final expression for the convective power -density transport 
is obtained: 

H k= c P T( ’ 7 ti » 

III. DERIVATION OF THE CONVECTIVE FLUX 
BY THE PHOTON TRANSPORT EQUATION 

In order to compute the power density transported by 
radiation, H rad , we apply the "gray spectrum" approximation 


5 
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with Rosseland's coefficient of opacity K . The source function 
is expressed as: 

B(t ) = -2-T 4 (t). 

To express the dependence of B on the optical depth 7, 
we use the equation of radiative transfer: 

4?r r 2 d(acT i ) 

^ = 3 Kp dr 

together with the geometrical relation: 

L(r) = 4jrr 2 H rad . 

The expression for optical depth is introduced: 

d7 = 7c pdr, 

and the Stefan-Boltzrrann constant: 

a = ac/4. 

These five relationships are combined: 

ac 6T* = 4 dT^ _ 4 dB 

H rad = “ 3 dr = " 3 dr = ' 3 d? ’ 

and the formula for hydrostatic pressure is used: 

dT = Kp Sd(lnP); 


therefore, 


H rad 


16 

3 



(2) 


The convective envelope is a fluid which possesses, in 
principle, all possible degrees of freedom. At large Reynolds 
numbers, very many of them come into play. The motion then 
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becomes so complicated that simple answers can be given to only 
those questions which concern the average behavior of the moving 
particles. Therefore, the equilibrium between the different ele- 
ments has to be described statistically. As in the kinetic theory 
of gases, the interest lies in the statistically most probable dis- 
tribution of the elements. 

Equation 2 for radiative power transport is derived from 
the equation of radiative transfer, which applies to all statistically 
controlled photon transport processes. Equation 1 for convective 
power transport is also based on statistical premises. The 
equation of radiative power transport will be used next as a base 
for deriving equations of convective power transport. 

The derivation of the equilibrium condition, which takes 
the place of Equation 1, is the same as the derivation of Equation 2 
except that convective parameters are substituted for the radiative 
parameters. The method of Schwarzschild is employed. 

The basic equation of radiative transfer, which is valid 
at every point within the Sun, has the form: 


91 

— cos 0 - 
3r 


91 

90 


-Si +7PI 


47 - 


Instead of considering the function I, which represents the 
distribution of the radiation in all directions from a point, we 
consider the first three moments of this distribution function: 


E(r) = i fldu 
H(r) = f Icos 0dw 
P rad (r) = c JlcoJddw 


radiant energy density 
radiation power flux-density 
radiation pressure. 
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Differential equations for these moments can be obtained 
by forming the corresponding mom nts of the equation of radiative 
transfer. The radiative transfer equation is multiplied by powers 
of cos B and integrated over all directions; then multiplying it 
by 1 and cos 0 , respectively, we obtain the first two moments: 


dH , 2 „ , _ _ . 

+ ~ H + ckp E - jp = 0 

^ + - <3P , -E) + H = 0. 

dr r rad c 


(3) 


These are two equations for the three radiation functions 
E, H, and P. To secure \ definite solution, we have to derive 
an additional relation between the moments. At a given point, the 
radiation field can be represented by the series: 

I = I 0 + I, cos 0 I ? cos 2 6 + 

As a result of the large absorption coefficient, this series 
converges rapidly in the solar interior. Therefore, the series 
will be restricted to the first two terms; it cannot be restricted 
to the first term only, because the radiation field would then be 
isotropic, without any net flux. If the first two terms of the series 
are introduced into the moment equations, the following results are 
obtained: 


E = 

4 7 r 
c 


H = 

4 i r 

i 

3 

i 


4 7T 

X o ' 

rad 

3 c 
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therefore, 


P rad= E ' 3 ‘ 


(4) 


Equations 3 and 4 form the required set of three equations 
for the three moments. 

The flux density H(r) may be represented in terms of L(r), 
the flux through the spherical area of radius r: 

H = L/47rr 2 . 


Sc far the derivation has been identical with that developed 
from the radiation equilibrium condition. However, instead of 
the Stefan -Boltzmann law being introduced, which is the next step 
in the radiative derivation, the analysis proceeds as follows: 

In order that the adiabatic equilibrium condition may be 
derived the negligible contribution of adiabatic mass expansion is 
recognized in defining the emission coefficient for the convecting 
mass: 

j = /Tp cc v T ; (c now denotes the 

velocity of sound). 

Substituting these relations for H and j into the first 
Equation 3 yields the simple form: 

E = c v pT . (5) 

Equations 4 and 5 yield the following expression for the 
radiation pressure: 

p rad = c v pT/3. (6) 

When Equations 5 and 6 are substituted into the second 
of Equations 3, the new expression for power flux-density transport 
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(convective) is obtained: 


H = 


cc v dT 
3? dr 


(7) 


When H is replaced by L/47rr , the basic equilibrium con- 
dition is obtained in the form: 

4iri* cc. 




dT 

dr 


(8) 


which is a new expression for the basic condition of convective 
equilibrium stated in the Introduction. 

According to Equation 5, the total internal energy of the 
convecting unit mass equals: 


U = c v T , 


which yields for the specific heat: 



Since Equations 1 and 7 are both valid expressions for con- 
vective power -density transport, they must be equal to each other. 
When the scale height S is set equal to the mixing length i in 
Equation 1 , the following relation is derived: 


v ( V - V') = 


2 c dT 

3 TRpT dr * 


IV. RESULTS AND CONCLUSIONS 

The heat flux in the convective layer of the Sun is derived 
by two different methods. The first method is conventional, and 
starts with the condition of thermodynamic equilibrium. The 
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second method is new, and is based upon the photon transport 
equation. The two formulas for the same quantity yield a relation 
between the temperature gradients occurring in the equations. 
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